Operator ordering and causality 
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It is shown that causality violations [M. de Haan, Physica 132A, 375, 397 (1985)], emerging 
when the conventional definition of the time-normal operator ordering [P.L.Kelley and W.H.Kleiner, 
Phys.Rev. 136, A316 (1964)] is taken outside the rotating wave approximation, disappear when the 
amended definition [LP. and S.S., Annals of Physics, 323, 1989 (2008)] of this ordering is used. 
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Introduction. — Causality in photocounting is a long- 
standing issue [1-;5| . No one really doubts that field radi- 
ation and detection is a causal process (cf, e.g., Refs. [f|), 
but to identify universal causal quantities measured by 
a macroscopic detector of arbitrary design remains an 
open question. Such quantities naturally emerge in re- 
sponse formulation of quantum electrodynamics (QED) 
[6j4l0|. In this letter we consider causality properties of 
these quantities and show that their use eliminates all 
causality issues from the theory. 

Referring the reader for details to the cited papers, 
here we only outline the key points. In Glauber's pho- 
todetection theory 11-^1, spectral properties of the 



detected field are determined by the quantum average 
(£(' (x, t)£( + > (x, t')\, where £^> (x, t) denote frequency- 
positive or negative parts of the Heisenberg field operator 
£{x,t) = £(+\x,t)+£<-^(x,t), and x comprises all field 
arguments except time. De Haan 0] and later Bykov and 
Tatarskii 0] remarked that (£ ^(x,t)£ W(x',tf)) is not 
a causal quantity, and suggested that the actual mea- 
sured quantity should be the time-normal (TN) average 



(T:£(x,t)£(x',t'):). 



(1) 



The symbol T- • ■ • '• denotes the TN operator ordering of 
Kelley and Kleiner (KK) [13, [T3|. Quantity © differs 
from 2Rc(£(~\x,t)£ ^ + ^(x',t')) by nonresonant terms. 
However, Tatarskii Q pointed out that quantity ([1]) may 
also exhibit acausal behaviour. Here we show that the 
idea of de Haan, Bykov and Tatarskii is basically correct. 
The problem is with the KK definition of the TN order- 
ing which needs generalisation beyond the rotating wave 
approximation (RWA). By using the amended definition 
of Ref. [7] all causality problems are eliminated. 

Characteristic properties of the TN ordering. — If 
quantities measured by a detector are universal, i.e., in- 
dependent of the detector, they should be identifiable in 
a QED theory of a field source. Let 



H.c. 



(2) 



be a Hermitian field operator, described by the standard 
bosonic creation and annihilation pairs aj_, a K , with u K (x) 



being complex mode functions. The field interacts with 
a quantum source according to the Hamitonian, 



t (t) = - J dxJ(x,t)E(x,t), 



(3) 



The nature of the mode index k, of the free current op- 
erator J(x, t) and of the source Hamiltonian (occuring 
implicitly) may be arbitrary. By definition, E(x,t) and 
J(x, t) are the interaction-picture (free) operators; the 
same operators in the Heisenberg picture are £ {x, t) and 
J(x,t). Free-field operators and the interaction Hamil- 
tonian is all one needs to construct the standard nonsta- 
tionary perturbation theory [pij . This makes the Heisen- 
berg operators formally defined, at least in perturbative 
terms. 

We define the (generalised) TN operator ordering pos- 
tulating that, 1) under the TN ordering, classical radi- 
ation laws apply directly to Heisenberg operators. For- 
mally, this is expressed by the relation between the TN 
averages of the field and current operators, 



(T:£{x 1 ,t 1 ) ■ ■ ■ £{x m ,t m ):) = J dx'idt'i ■ ■■dx' m dt' m 

x A R (a;i,x' 1 ,ii - t' x ) ■ ■ ■ A R (x m ,x' m ,t m - t' m ) 

x (T:J(x' 1 ,t' 1 )---J(x' 1> t' m ):), (4) 

where Ar (x, x' , t — t') is Kubo's linear response function 
of the free field, 



A R (s, x',t-t') = -6{t - t') [E(x, t), E(x', t')] 



= -26(t-t')Im^2u K (x)u* K (x')e 



-iuj K (t— t ) 



(5) 



The averaging in (U) is over the Heisenberg (initial) state 
of the system. The latter is assumed to factorise into 
the vacuum state of all oscillators (denoted |0)) and an 
arbitrary state of the quantum source. We retain the con- 
ventional notation for the TN ordering, implying that the 
KK definition is a resonance approximation to it (which 
is indeed the case, see below). 



2 



Linear response functions of free bosonic fields and of 
the corresponding classical ones are identical [6j. Hence 
eq. Q exactly emulates a classical formula for stochas- 
tic averages of a field radiated by a random current into 
empty space (vacuum). Such formula is found from eq. 
(U) by dropping hats and replacing TN averages by clas- 
sical stochastic averages. Note that condition 1 warrants 
explicit relativistic causality in radiation and propagation 
of the field but says nothing about causality properties 
of the TN ordering itself. 

Furthermore, 2) for free electromagnetic operators, the 
TN ordering coincides with the conventional normal or- 
dering [l3j ]. This is a consistency requirement: with 
J{x,t) = eq. (H|) turns into, 



(0\T:E(xi,h)---E(x m ,t m y.\0) = 0. 



(6) 



This relation is enforced by condition 2. 

Time-normal ordering beyond the RWA. — Character- 
istic conditions 1 and 2 are supplemented by formal ones: 
3) the operation T: ■ ■ ■ : is polylinear, 4) any commuting 
factor may be taken out of the T: • • • : symbol, and 5) 
T:l: = 1. The TN product of m arbitrary operators 
Ak(t), k = 1, ••• , m, obeying conditions 1-5 [22[, is de- 
fined as @,H, 

T :Ai{h) ■ ■ ■ Am{t m ): 

/m 
dt' 1 .--dt' m T c Y[[6^(t k -t' k )A k+ (t' k ) 
k=i 

+ 5(-\t k -t' k )Ak-(t' k )]. (7) 

The kernels are the frequency-positive and negative 
parts of the delta- function, (cf. Ref. Q , appendix 2) 



5^(t) = ± 



27n(i=R0+)' 



(8) 



and the Tc, or closed-time-loop, ordering [15| is a way 
of writing the double-time- ordered operator structure 
(known, e.g., from the photodetection theory 12|, [l3j), 



TAxih) ■ ■ ■ A m (t m ) T&(ti) ■■■B n {t' n ) 

= Tcii-(ti) ■ •■i m -(t m )Bi+(ti) • --Bn+it'J, (9) 

where T is the standard time ordering of operators and 
T is the "reverse" ordering. The ± indices serve only for 
ordering purposes and otherwise should be disregarded. 

Polylinearity of (|7|) (property 3) is inherited from the 
Tc-ordering. Property 4 is a consequence of the relation, 



+ <5(->(t) = 6(t). 



(10) 



Property 5 is a specification of ([7|) for m = 0. Property 
2 is verified in Ref. Q . The one most difficult to demon- 
strate is property 1. Its formal proof 0, [l(| employs 
heavy-duty machinery of quantum field theory. However, 



the basic idea is fairly simple. The starting point of the 
proof is the wave quantisation formula 



17] 



[E(x,t),E(x',t')] 



-ih[A R (x,x',t - t') 
- A R (x',x,t' 



t)]. (11) 

It is found inverting Kubo's formula ([S]). Equation (fTTj) 
induces restructuring firstly of Wick's theorem, and con- 
sequently of the whole standard perturbative approach 
of the quantum field theory. Equation (|4]) is a rigorous 
form of eq. (74) in paper 

To extend to fermionic operators one may assume 
that they always occur multiplied by generators of an 
auxiliary Grassmann algebra 

[Emm, 



A k (t k ) = 7kFk{tk)- 



(12) 



Such combinations behave under orderings as bosonic 
operators Q- Quantum fields are included by making 
all field "labels" explicit. For instance, in spinor elec- 
trodynamics, A k {tk) = ^ fc (?'fc,tfc),7fc^o' fc ( T 'fc,ifc), where 
Afj_(r,t) is the 4- vector electromagnetic potential and 
ip a (r,t) is the electron-positron field. Other quantum 
fields may be included similarly. Equation ([7]) thus ex- 
tends the KK definition in two ways: beyond the RWA 
and to arbitrary quantised fields including fermions. 

Time-normal ordering under the RWA. — Separation 
of the frequency-positive and negative parts of a func- 
tion, occuring in Kelley-Kleiner's definition of the TN 
ordering, may be written as an integral transformation, 
(see, e.g., Ref. [7[, appendix 2) 



f (± \t)= dt'8^\t-t')f{t'). 



(13) 



In ([7]), the Tc-ordering applies to entire operators, and 
not to their frequency-positive and negative parts (i.e., 
Tc first, W second). The KK definition emerges by 
changing the order of operations first, Tc second) 
Q. This results in a resonance approximation to ([7]): 



T -.Aiih) ■ ■ ■ Am(t m ): 



Tcl[[A k X\t k ) + Ai-J(t k )], (14) 



fc=i 



which indeed coincides with the KK definition. 

Parametric oscillator. — The simplest example of a 
system which exhibits a causality violation with defini- 
tion (fT4|) is the parametric oscillator, 



H(t) 



p 2 (t) muj 2 {t)x 2 (t) 
2m + 2 



(15) 



Its frequency loq is reduced by half at t = and restored 
to its initial value at t — 2Tq = 47r/wo: 



= ( IP ° <t<2To < 
I w , t < or t > 2T . 



(16) 
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FIG. 1: The time-normal average (0\T:p(t)p(t):\0) (solid 
line), the KK approximation to it (dashed line) and the aver- 
age (0|p 2 (f)|0) (dash-dotted line) in units of fimwo- The graph 
of (0\p 2 (t)\0) is shifted vertically by Smwo/4. 



The initial (Heisenberg) state of the oscillator is vacuum 
(defined with respect to ojq). For p(t) we have, 



Pit) 



u t mw i U) t 
pcos — — sm — > < i < 2T , 

pcosu>Qt — mjjjQX sinwoi, t < or t > 2Tq. 



(17) 



The quantity (fS\T-p{t)p{t'): |o) for t — t' calculated using 
eqs. (J7J) and (|T3|) is drawn in Fig. Q] with solid and dashed 
lines, respectively. The former is zero for t < 0, whereas 
the latter is nonzero for all times. Using eq. (fT4)) indeed 
resulted in a causality violation. 

Used here solely as an illustration, the parametric oscil- 
lator appears to be of much interest by itself. The physi- 
cal motivation for calculating (p\T:p(t)p(t'): |0) is that it 
describes a continuous measurement of fluctuations of the 
quantised momentum by means of electromagnetic inter- 
action. The quantum current associated with the particle 
is J(t) — ep(t)/m, where e is its charge. Radiation of 
this current is observed by a remote detector. In essense, 
Hamiltonian (TT5j) is postulated as an effective one for the 
Heisenberg current operator J(t) = ep(t)/m. Figure Q] 
demonstrates two things. Firstly, the time- normal or- 
dering matters. For comparison, we also draw in Fig. Q] 
the quantity (0|p 2 (i)|0) (dash-dotted line) @. Clearly 
the exact and the KK results for (0|T:p(t)p(t'):|0) are 
much closer to each other than to (0|p 2 (i)|0). Secondly, 
all differences are limited to short transients at t = 
and t = 2T . Up to a vertical shift, the beyond-the-RWA 
(0\T:p(t)p(t'):\0) and (0|p 2 (i)|0) quickly settle into the 
same pattern; for Tq/2 <t < 2Tq they are indistinguish- 
able. Furthermore, the reader may be surprised — and 
even disturbed — by the transient seen in Fig.[T]for times 
beyond 2Tq. It should not be overlooked that the oscil- 
lator returns to its vacuum state at t = 2T Q ■ As sug- 
gested by A.Kaplan 11911 ■ this transient may be a toy case 
of the Unruh effect (2Cj. Indeed, since (0\J{t)\0) = 0, ra- 



diation of the current j(t) is spontaneous emission. The 
transient reflects its non-instantaneity, which in turn is 
due to the time-energy uncertainty relation. However, to 
claim results in Fig. [1] as physical, one has to account for 
the radiation friction and its fluctuations disregarded in 
the effective Hamiltonian. We return to this discussion 
elsewhere. 

"No-peep-into-the- future" theorem. — We now show 
that a TN product depends on the operators it comprises 
only for times not later than its latest time argument. 
The question of causality of TN products thus reduces to 
that for quantum equations of motion. 

Proof of this theorem is a simplified version of the proof 
of causality of quantum response functions in Ref. Q- 
We break the integration in ([7]) into m domains, labelled 
by k = 1, • • • , m, such that, 



t> ... t> t' 



(18) 



Within the fcth domain, the &th factor in the product in 
([7]) is transformed as follows, (cf. eq. (QI 



dt' k [^-\t k -t' k )A-{t' k ) + 5^\t k -t' k )A + {t' k ^ 

dt' k A + {t' k )[5(-Ht k -t' k ) + 6^(t k -t' k )] 

= A+(t k ). (19) 



Indeed, the latest operator in eq. (|9]) is positioned be- 
tween the T and T-ordered products irrespective of its 
C-contour index (cf. eq. (110) in Ref. (tJ). The integra- 
tion over remaining variables is restricted to, 



(20) 



making the "no-peep-into-the-future" theorem evident. 

Relativistic causality. — As a comparatively simple ex- 
ample, consider the pair time-normal product, 



T :Ai{r 1 ,t 1 )A 2 {r2,t 2 ) 



dt[ 



dt' T, 



2 1 C 



x - ii)ii+(ri,ti) + - 4)A-(ri,4)] 

x [*(+> (t 2 - t' 2 )A 2+ (r 2 ,t' 2 ) + *(-> (t 2 - t' 2 )A 2 - (r 2 , t' 2 )} . 

(21) 

Hereinafter we assume that t\ > t 2 . In the spirit of 
the above proof, we break the integration into the two 
domains, t\ > t' 2 and t' 2 > t\. In both domains, time 
ordering is performed explicitly. Rearranging the terms 
and using (|10p we find, 



T :A 1 {r 1 ,t 1 )A 2 (r 2 ,t 2 ): 

= t dt'^+Xh-t'JA^^MriA) 

J — OO 

+ 6^ (t 2 -t' 2 )A 2 (r 2 ,t' 2 )A 1 (r 1 ,t 1 )]+ {1^2}. (22) 
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verified for a time-normal average of two operators, while 
extention to more than two operators and, much more 
importantly, to renormalised theories remains subject to 
further work. 
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FIG. 2: Space-time geometry of the integral in eq. (|21|l . 
Points r\ , t\ and r 2 , t 2 are shown as bold dots, and their past 
light cones as slanted lines. The thick vertical line represents 
the "integration path" r 2 ,t' 2 , —00 < t' 2 < t\. It consists of 
two sections, positioned inside (—00 < t' 2 < t, solid) and out- 
side (i < t' 2 < ti, dashed) of the past light cone of ri,ti. 
Other lines guide the eye. 

The integral here is illustrated in Fig. [21 where the space- 
time points ri,t\ and r2,^2 are drawn as bold dots, and 
the domain ri,t[, —00 < t\ < t 2 ("integration path") as 
a thick vertical line with solid and dashed sections. 

Now, what kind of "no-peep-into-the-future" theorem 
would one expect in relativity? Assuming that the depen- 
dence of operators on various perturbations is relativis- 
tically causal, it is sufficient to assume that space-time 
arguments of the operators which a time-normal product 
comprises are confined to the union of the past light cones 
of its arguments. Without assumptions about quantum 
dynamics, eq. (|22|) explicitly violates this condition. So, 
in Fig. [21 the integration path extends into the future be- 
yond the point r 2l t 2 . The minimal dynamical assump- 
tion one has to make is that two operators commute if 
their arguments are separated by a space-like interval (it 
is one of Wightman's axioms |21|). Such commutativ- 
ity holds for A\{r\,t\) and A2(r2,t' 2 ) with t < t' 2 < t\, 
where t is the time when the integration path exits the 
past light cone of r\,t\ (Fig. [2]). With this observation, 
the contribution from the dashed section of the integra- 
tion path in Fig. [2] reduces to, (cf. eq. (fTO]) ) 

= 9(t 2 -f)i 2 (r 2 ,*2)ii(ri,ti). (23) 

The offending contribution from t' 2 > t 2 cancelled. Simi- 
lar arguments apply to the second term in (|22l) [25j |. 

Conclusion. — We have demonstrated, both by exam- 
ple and by a formal proof, that causality violations, oc- 
curing if Kelley-Kleiner's definition of the time-normal 
operator ordering is taken outside the rotating wave ap- 
proximation, are eliminated by putting the amended def- 
inition of Refs. 0, H| to use. Relativistic causality was 
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